Contact induced spin relaxation in Hanle spin precession measurements 
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In the field of spintronics the "conductivity mismatch" problem remains an important issue. Here 
the difference between the resistance of ferromagnetic electrodes and a (high resistive) transport 
channel causes injected spins to be backscattered into the leads and to lose their spin information. 
We study the effect of the resulting contact induced spin relaxation on spin transport, in particular 
on non-local Hanle precession measurements. As the Hanle line shape is modified by the contact 
induced effects, the fits to Hanle curves can result in incorrectly determined spin transport properties 
of the transport channel. We quantify this effect that mimics a decrease of the spin relaxation time 
of the channel reaching more than 4 orders of magnitude and a minor increase of the diffusion 
coefficient by less than a factor of 2. Then we compare the results to spin transport measurements 
on graphene from the literature. We further point out guidelines for a Hanle precession fitting 
procedure that allows to reliably extract spin transport properties from measurements. 
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I. INTRODUCTION 

New concepts like the spin transfer torque, the trans- 
port of spin information over long distances and the 
prospect of spin field effect transistors keep spintronics 
an inspiring fielcP^l. But before new types of spintronic 
devices can be build we need both materials that 
efficiently generate spin currents as injector and detector 
electrodes and materials with long spin relaxation 
lengths (A) and times (r) to transport the spins with 
only little losses. 

While ferromagnetic metals can spin polarize currents 
and are therefore used to inject and detect spins, semi- 
conductors offer low spin relaxation which make them 
good candidates to be used as transport channels. One 
of the main challenges when combining the two t ype s of 
materials is the "conductivity mismatch" problem^. As 
the electrical resistance in the ferromagnetic electrodes 
is in general lower than in the semiconducting transport 
channel, the injected spins tend to be reabsorbed by the 
leads and loose their spin orientation. 
Graphene being an intermediate between metal and 
semiconductor systems is a prototype example for the 
conductivity mismatch, as graphene based devices can 
be well described following simple spin diffusion theory. 
Here, the long spin relaxation lengths of several fim 
measured at room temperature are already promising 5 , 
but still stay behind the theoretical prospects based 
on the high mobilities combined with weak spin orbit 
coupling and low hyperfine interactions^. While some 
research aims to understand the spin relaxation mech- 
anism in graphene^^ and to understand the influence 
of the direct environment of the graphene transport 
channeP^H^ the conductivity mismatch can play an 
important role in the origin of spin relaxation in the 
measured devices. To prevent this mismatch high 
resistive barriers be tween the contacts and the graphene 
channel are includecPE^HH] 



The most common and reliable way to probe spin 
transport properties is by perfo rming measurements in 
the non-local spin- valve geometry * 2 * 5 * 24 ! because it enables 
to separate spin and charge currents, avoiding spurious 
effects^. Popinciuc et alP^D describe, in agreement 
with Takahashi and MaekawaP^, that the measured 
amplitude of the spin signal in the non-local geometry 
is strongly reduced for low contact resistances RcW^. 
To quantify the effect Popinciuc et al. introduce the so 
called i?-paramctcr that is defined for a 2-dimensional 
channel by R = [Re / R sq )W , where R sq is the s quare 
resistance and W the width of the diffusive channel!- 19 * 28 ! 
In this article we start by summarizing the dependence 
of the non-local amplitude on the contact resistance 
discussed in Ref. [19] and then we are going to focus on 
determining how the Hanle spin precession is influenced 
by low contact resistances. We discuss that not only 
the amplitude but also the shape of Hanle precession 
curves is changed for low values of the i?-parameter 
(corresponding to low contact resistances) and simulate 
Hanle measurements including contact induced relax- 
ation. We quantify the contacts' influence by fitting 
the data with the standard Hanle formula without 
taking contact induced effects into account, assuming 
R — » oo. Note that fitting with the standard Hanle 
formula is the common method to analyze experimental 
spin precession data in almost all published work. The 
difference of the extracted spin relaxation time T^ lt and 
diffusion coefficient D^ %t to the parameters used in the 
simulations is quantified and we compare these results 
to data obtained on graphene spin- valve devices where a 
reduction of t was reported for low contact resistances^. 
Finally, we point out how to extract correctly the spin 
transport properties from Hanle precession measure- 
ments by excluding spurious background effects. 
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FIG. 1: (Color online) Sketches of spin diffusion through a 
diffusive channel with a spin injector and detector separated 
by a distance L in non-local geometry for (a) high and (b) low 
contact resistances. The width of the contacts in x-direction, 
Lc, is small compared to L or the spin relaxation length A. (c) 
The spin chemical potential fis indicating the spin accumu- 
lation below the injector electrode and the exponential decay 
of the spin signal (red dotted curve). The spin accumulation 
influenced by the contact induced spin relaxation is denoted 
by the black solid curve. 



II. CONTACT INDUCED SPIN RELAXATION 

Fig. [I] (a) presents a sketch of the non-local measure- 
ment geometry with an injecting electrode at x = L on 
the right side, sending a charge current via a resistive 
barrier^ into the channel to the right side end and a de- 
tecting electrode on the left side (at x = 0), measuring 
the voltage difference between the contact and the left 
side end. As the electrodes are ferromagnetic, the inject- 
ing electrode induces a spin imbalance that accumulates 
below the electrode and diffuses away from it in both pos- 
itive and negative x-direction of the channel (red dotted 
curve in Fig. [I](c)). The second ferromagnetic electrode 
detects the spins at x — and the measured voltage is 
normalized with the injected curre nt to ob tain the non- 
local resistance R n i that is given by * 19 * 26 * 2 ^ 



±- 



(2i?/A) 2 exp(-L/A) 



2W (1 + 2i?/A) 2 - exp(-2L/A) 



(1) 



The model leading to this result is based on the one- 
dimensional description of a diffusive channel with an 
injector and detector on distance L with P the polariza- 
tion of the contacts and A = V Dt the spin relaxation 
length in the channel with the diffusion coefficient D. 
The width of the contacts (Lc) is considered to be negli- 
gible compared to L and .)pS. Also we assume 1 — P 2 « 1 
(applicable to graphene devices where P < 309cP^) and 



are considering an infinite homogeneous transport chan- 
nel. The effect of an inhomogeneous transport channel 
is discussed elsewhere^. The i?-parameter is being cal- 
culated using the contact resistance of the injector and 
detector. In case R is not the same for the two electrodes 
an effective single i?-parameter can be calculated using 
1/R e ff ~ (l/i?i + l/i?2)/2 with the i?-parameters of the 
injector and detector R\ and i?2 (see Appendix |A|). The 
meaning of the i?-parameter gets clear when it is normal- 
ized with the spin relaxation length A. The normalized 
value corresponds to the ratio of the contact resistance 
and the spin resistance of the channel R s ch = R sq X/W 
so R/X = Rc/R s ch - Hence, R/X describes the ratio of 
spins diffusing through the channel and relaxing, versus 
those being reabsorbed by the contact, making it a good 
measure for the influence of the contacts. 
Eq. |T]) shows that the spin signal R n i has a maximum 
for high contact resistances (R — > oo) and is reduced 
for low R- values. A significant change is observed for 
R/X < 1 (Fig. [2] (b), inset). On the other hand the am- 
plitude of the signal is reduced with increasing L from 
a maximum at L = 0. The characteristic length ratio 
of the system is L/X. While the effect on the normal- 
ized amplitude (R n i /RnT^°° W1 th the amplitude without 
contact induced effects R^ 00 ) is smaller for short dis- 
tances between injector and detector electrode it stays 
approximately constant for L/X > 1. Popinciuc et al. 
discuss in detail the effect of low contact resistances on 
the measured non-local amplitude but, while included in 
the model, the effect on the Hanle curve is only discussed 
qualitatively^. In the following we are going to present 
a quantitative analysis of the influence of low contact re- 
sistances on Hanle measurements. We show that the ex- 
tracted spin transport properties of the transport channel 
can be limited by the contact induced relaxation and are 
therefore incorrectly determined when low R measure- 
ments are analyzed without considering the influence of 
the contacts. 

Fig. 2] (a) shows Hanle precession data that was simu- 
lated for different values of R/X with L/X — 1 using the 
model system of Fig. [T] (a) described in Ref. '19, Note 
that the amplitude of the Hanle curves is normalized at 
B = 0, which is necessary as the amplitude scales with 
(R/X) 2 for R/X <C 1 and changes by 5 orders of magni- 
tude between R/X = 0.001 and R/X = 10. A significant 
change in the Hanle shape is visible in Fig. [2] (a), point- 
ing to an effective change of the spin transport proper- 
ties. The strongest change in the shape is seen between 
R/X = 0.01 and R/X = 1 while the curve shape is sat- 
urating for both small and large R/X values denoting 
spin transport limited by the contacts or by the prop- 
erties of the channel, respectively. Fig. [2] (b) shows a 
similar dataset for L/X = 10. We also see a change in 
the Hanle shape but the effect is much weaker for this 
larger distance of the injector and detector. Remarkably, 
in both cases the curves stay in the characteristic Hanle- 
like shape for all R/X. Therefore it is possible to fit the 
data using the solutions to the Bloch equation d 29 * 30 1 that 
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FIG. 2: (Color online) Simulated spin precession curves for 
different values of R/X with (a) L/A = 1 and (b) L/A = 10. 
For the simulations we use D = 0.01 m 2 /s, r = 100 ps, W = 
1 )im and R sq = 1 kfi (representative of graphene devices) 
with contact resistances between 1 and 10 4 fl. The amplitude 
of the curves is normalized for clarity with R n i(B = 0). The 
inset in panel (b) shows R„i from Eq. as a function of R/X 
for L/A = 1, normalized by R n i(R — > oo) (black solid line) 
and the asymptote oc R 2 (red dashed line). 



do not take the effect of the low resistive contacts into 
account (see Appendix |Tj| 



dt 



DV 2 fT s 



Ms 



+ bT L x (f s . 



(2) 



Here ail is the Larmor frequency oil = g\iB /h B, with the 
gyromagnetic factor g (g- factor, g ps 2 for free electrons), 
the Bohr magneton fig and the magnetic field B. By 
fitting simulated data without taking the effects of the 
contacts into account, we can determine what happens 
when one fits the data obtained in samples with corre- 
sponding R- and L-values in the standard manner.^ 
The results from these fits are presented in Fig. [3] Note 
that while the simulations were performed with D = 
0.01 m 2 /s, r = 100 ps they do not depend on the spe- 
cific value of D and r. Hence, we get the same results 
for different D and r resulting in the same A = \J Dt 
as the fitting results depend only on the ratios R/X and 
L/A. The graphs show the fitted values r^ 1 * and D* %t 
normalized by the actual values for the channel r and 
D as a function of R/X for different L/A (Fig. [3] (a) and 
(c)) and as a function of L/A with different R/X (Fig. [3] 
(b) and (d)). While all values converge for high R/X to 
the intrinsic values, we see a strong decrease of r^ lt and 
a moderate increase of D^ %t for small R/X. Looking at 
Fig. [3] (a) in more detail we see that the decrease in T? lt is 
the strongest the shorter the distance L between injector 
and detector relative to A. We also see that the values 
saturate for small values of R/X as already perceivable 
in Fig. [2] (a). In this limit the effect of the contacts is 
maximized, t' %t shows changes of nearly up to five orders 
of magnitude which means that in a measurement with 
parameters of R/X = 0.001 and L/A = 0.01 we would 
underestimate r by a factor of 5 x 10 4 . 
The length dependence of the effect is more clearly pre- 
sented in Fig. ^ (b) where the r^ 1 * data is plotted as a 
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FIG. 3: (Color online) The change in r fit and D fit fitted 
for different L/A as a function of R/X ((a) and (c)) and for 
different R/X as a function of L/A ((b) and (d)). For small 
values of L/A the fits become insensitive to the specific value 
of the diffusion coefficient, resulting in the non-smooth be- 
havior shown for L/A = 0.1 in panel (c). Therefore the data 
for L/A < 0.1 is not shown in panel (c) and (d). 



function of L/A for different R/X. Here we see that while 
the decrease of r^ lt is stronger for shorter distances the 
effect gets negligible for L/A > 10. That means that 
contact induced effects can be circumvented by measur- 
ing on a longer distance. This is only limited by the 
reduced measured amplitude for longer distances L (see 

Eq. m. 

Fig.klnc) and (d) show the same plots for D^ lt . Also here 
we see the strongest effect for small R/X and L/A and no 
significant change for R/X = 100 or L/A = 100. On 
the other hand, the values for D' lt show a much weaker 
change than the values of T^ lt and the change is directed 
in the opposite direction than the change of r^ lt . Similar 
to t-^*, the D-f lt values also seem to saturate for small 
R/X and the changes are less than a factor of 2. 
While most curves presented in Fig. [3] have a smooth 
shape and a continuous change with L/A and R/X, the 
data for D^ lt shows for values of L/A < 0.1 combined 
with values of R/X < 1 a non-smooth behavior. This is 
related to the fact that the diffusion in the channel gets 
for small L dominated by the contact induced effects for 
short distances and low contact resistances and the shape 
of the Hanle curves gets strongly influenced. The spin ac- 
cumulation has no significant decay between the injector 
and the detector electrode so the system becomes sim- 
ilar to 3-terminal Hanle precession^. As a result, the 
fits become insensitive to the specific value of D and one 
can only determine 7^ (see Appendix [bJ . Therefore we 
omitted the data for D flt for L/A < 0.1 in Fig. [31(c) and 
(d). 
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Note that in the limit L/X <C 1 and R/X < 10 the values 
for r^ lt saturate as they are dominated by the contact 
induced effects and can be described by a basic formula 
related to the back diffusion of the spins into the contact 
(see Appendix [C]) . 



III. DISCUSSION 

Fig. [3] shows clear trends for t^ %1 /t and D^ lt jD as a 
function of R/X and L/X. We are going to discuss in 
the following how to understand the physics behind the 
presented results. The sketch in Fig. [I] (a) presents the 
spin injection and detection for high contact resistances, 
e.g. due to tunnel barriers between the channel and the 
contacts. Here the spin diffusion in the channel remains 
undisturbed and the injected spins diffuse freely through 
the channel before being detected by the spin sensitive 
detector. In this way, measurements detect the intrinsic 
spin transport properties of the channel and the simple 
exponential decay of the spin signal (red dotted curve in 
Fig. [l](c)) is obtained. 

In the case of low contact resistances the spin transport 
is influenced both at the injector and at the detector 
electrodes. When diffusing through the channel the low 
resistive detector has a high probability of detecting the 
spins as soon as they are near the contact as it acts as 
a spin sink (Fig. [I] (b)). Therefore the effective traveling 
time is reduced and the measured diffusion coefficient en- 
hanced as D — L 2 /2t]j where Tp is the diffusion time for 
the length L. At the same time the proximity to the low 
resistive contacts also causes spins to relax, which reduces 
the relaxation time. The extra relaxation is depicted by 
the kink at the detector in the black solid curve in Fig. [T] 
(c), describing the decay of the in general reduced spin 
accumulation in the system. 

Fig.|3](b) and (d) show a reduction of the contact induced 
effects for larger L/X. This can be easily understood by 
the fact that for a longer distance between the electrodes 
the ratio of the time the spins stay in the channel com- 
pared to in close proximity to the contacts grows result- 
ing in relatively less influence of the contacts on the spin 
transport. 

Fig. H (a) shows the spin relaxation length X^ lt resulting 
from the fitting results presented in Fig. [3] for L/X > 0.1. 
The shape of the A^*-curve is comparable to the behav- 
ior of T^ %t . This is due to the fact that X$ lt is mainly 
influenced by the spin relaxation time t' lt with a change 
of up to a factor 1000 (for L/X > 0.1). As D* u shows 
only a change of less than a factor 2 we get a maximum 
reduction of X' lt by a factor 25. 

This X^ lt value would be used in the analysis of a mea- 
surement to calculate the polarization P. If we take 
the amplitude simulated with Eq. |l]) (inset, Fig. [2](b)) 
and assume spin transport without contact induced spin 
relaxation (Eq. Q for R — > 00) we extract the effec- 
tive polarization P?" with R n i{X flt , P flt , R -> 00) = 
R n i(X, P, R). The resulting value is up to 500 times re- 
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FIG. 4: (Color online) (a) The change in X flt /X calculated 
using T^ lt and D^ lt from Fig. [3] (b) The effective polarization 
P flt normalized with the actual polarization P. The values 
are plotted as a function of R/X for different L/X. The values 
for L/X = 100 and 10 overlap in both panels. 



duced for small values of R/X compared to the real P 
value (Fig. [4] (t>) ) . Note that the largest change in Pf %t 
compared with P is observed for long distances. 
After discussing the effects observed in the simulations 
let us have a look on measurements on real devices us- 
ing graphene as the transport channel. In spin transport 
samples in graphene it is difficult to produce high resis- 
tive contacts and to control the quality of the contact- 
graphene interface. So a data set with similar quality 
samples with only a change of the contact resistance 
is difficult to produce. On the other hand in a single 
device the quality of the contacts is most of the time 
comparable. Therefore it is relatively easy to check the 
length dependence of the spin transport properties in this 
kind of system, assuming similar i?-parameters for all 
electrodes. Two sets of data obtained on two different 
graphene devices with three different injector-detector 
distances are presented in the work by Wojtaszek et al. 
in Fig. 4 (a) and Fig. 5 (a) of the supplementary infor- 
mation of Ref. |33l In both cases a minor increase of D 
is reported when measuring on a shorter distance and 
in the first case also a minor decrease of r, pointing to 
weak but apparent contact induced relaxation^. With 
R > 3 fim and A ~ 5 /jm the measurements were also 
performed in a regime where one would expect this kind 
of weak contact induced effects as L/X « R/X s» 1 (see 
Fig. |p. 

Han et al. present a study of the dependence of the 
spin transport properties on the quality of the resistive 
barrier between the graphene channel and the contacts 
in Ref. 1201 They show that between tunneling injection 
of spins and the injection with transparent contacts the 
measured spin relaxation time decreases while the diffu- 
sion coefficient is increased in agreement with our sim- 
ulations' results. On the other hand the results for a 
"pinhole" barrier with intermediate resistance present an 
intermediate spin relaxation time but also a reduced dif- 
fusion coefficient. While the spin relaxation time fits into 
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the expectations for an intermediate contact resistance, 
the reduced diffusion coefficient cannot be explained by 
the contact resistance but has to be related to a lower 
quality sample or other effects. 

Our model also points to the fact that the recent re- 
ported differences between the results for the spin relax- 
ation length, based on the analysis of 4-terminal non- 
local Hanle precession measurement^] and based on the 
analysis of the magnitude of spin- valve measurements in 
local 2-terminal geometry with very high contact resis- 
tances (Rc > 1 Mf2)p3 cannot be explained by contact in- 
duced relaxation. If one would measure with the configu- 
ration of Han et alP^with L — 5.5 /im and R ~ 200 /im a 
material with a spin relaxation length of ~ f 00 /im and a 
spin relaxation time of ~ 100 ns as reported in Ref.l2"2"lone 
would only see a reduction of the fitted spin relaxation 
time by a factor of t^ u /t w 1/3 (see Fig. (a)) lead- 
ing to a reduced spin relaxation length of X^*/A ~ 1/2 
(see Fig. [4] (a)) as one would have L/X « 0.05 and 
R/X sa 2. Therefore the standard Hanle analysis would 
yield X^ lt ps 50 /im and t* u ~ 30 ns but Han et al. report 
X fit « 2.5 /im and r fit « 0.5 nsP. With X flt w 2.5 /mi 
Han et al. are in the regime of negligible contact induced 
relaxation with L/X ~ 2 and R/X ~ 80, so the differ- 
ence in the measured A is not based on contact induced 
relaxation but has to be related to other effects. Even 
for a spin relaxation length of A = 20 /im it would be 
L/X ~ 0.25 and R/X s» 10 for the system of Han et al. 
and they would be able to measure this A without sig- 
nificant influence of the contacts (see Fig. [3] (a)). Such 
strong differences of the spin signal magnitude between 
non-local and local configuration as between Refs .12171 and 
[22] have also been observed in traditional semiconductors 
like silicon in the non-locaP^ and 3-terminaP^ configura- 
tion. 



IV. GUIDELINES FOR A GOOD AND 
RELIABLE HANLE FIT 

In this paper we discuss how Hanle measurements are 
influenced by contact induced relaxation that can lead 
to incorrectly determined spin transport properties of 
the channel. Independently from that, the fitting proce- 
dure can also give incorrect results for the spin transport 
properties when performed incorrectly. In this section we 
are therefore commenting on typical pitfalls in analyzing 
Hanle precession data. 

The fit to a Hanle curve is unambiguous if performed in 
the right way. Fig. [5] (a) illustrates how a fit can still 
give wrong results on the example of a Hanle precession 
fit when assuming a wrong background resistance. The 
background resistance is represented by the R n i(B — > oo) 
value and is the fitting baseline. Fig. [5] (a) shows a fit 
to the central peak of a Hanle curve (without contact 
induced effects) with a baseline shifted by +5% of the 
amplitude. The fit results in an increase of r by > 10% 
and of D by > 45% and therefore a misestimation of A by 
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FIG. 5: (Color online) (a) The influence of a baseline shift 
shown by means of a Hanle spin precession curve, shifted +5% 
of the precession amplitude upwards (black solid curve) and a 
fit assuming no shift (red dotted curve). The baseline of the 
fit is therefore at R n i = while the baseline of the data is de- 
noted by the black dashed line. The same Hanle curves on a 
larger B-field range are shown in the inset. A clear difference 
is visible for |B| > 0.3 T. (b) The change of the diffusion co- 
efficient and the spin relaxation time resulting from data with 
a baseline shift and fits assuming no baseline shift. The pre- 
sented data was simulated using D — 0.01 m 2 /s, r = 100 ps 
and L = 1 jjxa.. 



more than 25% compared to the values used to simulate 
the data. However, when fitting the curve with these val- 
ues the fit presents itself faulty when including the high 
field tails of the curve as shown in the inset of Fig. [5] 
(a). Fitting to high B-ficld values gives therefore a good 
indication of the quality of the fit. However, this identifi- 
cation of a bad fit can be partly masked by data noise in 
combination with anisotropic magneto resistance effects 
or the out-of-plane tilting of the magnetization of the fer- 
romagnetic electrodes at high field values, adding an ad- 
ditional background resistance^. Another indication of a 
good fit is the fitted curve reproducing the "shoulders" of 
the measured curve, where R n i has a minimum (for par- 
allel alignment of the injector and detector electrode). 
This is obviously not accomplised in the presented case 
(Fig. [5] (a), inset). The larger the ratio L/X the more 
pronounced are the shoulders, so measuring on a longer 
distance enhances the reliability of the fit . 
While measuring to high magnetic field values to deter- 
mine the background resistance is in any case advisable, 
there is a way to avoid such spurious background effects 
in a fit. Measuring the spin precession both for parallel 
and for antiparallel orientation of the electrodes, and sub- 
tracting the signals from each other removes most spu- 
rious (not spin rel ated) bac kground effects as done in 
several recent work^ESMET] gy taking the mean of the 
parallel and the antiparallel measurements, one can also 
extract the -B-field dependent background resistance. Fi- 
nally, a minor error in a fit can also occur if the magnetic 
field values are not properly calibrated. The effect of a 
correction factor for the magnetic field value is the same 
as the effect of a changed g- factor as u>i oc B and wj, a g. 
A wrong .B-field calibration is therefore linearly passed 
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Appendix A: The i?-parameter for dissimilar 
contacts 



We discuss the effect of low resistance contact induced 
spin relaxation on Hanle precession data and quantify the 
misinterpretation of spin transport properties in a trans- 
port channel that can arise from this effect. As fitting 
Hanle curves is a common way to extract spin trans- 
port properties we use the model presented in Ref. [19] to 
simulate Hanle measurements and fit the data using the 
standard formula, neglecting the contact induced effects. 
The observed rescaling of the spin relaxation time and 
the diffusion coefficient only depend on the ratios R/X 
and L/X and the fitting results show that a strongly de- 
creased t' lt by up to nearly five orders of magnitude and 
a moderately increased D' zt by less than a factor of 2 can 
be observed for small R/X and L/X. On the other hand 
large values for both R/X or L/X show a convergence of 
T^ lt and D? lt on the undisturbed values r and D, inde- 
pendent on L or R, respectively. This shows that the spin 
relaxation induced by the contacts can in principle be 
avoided when measuring on a longer distance. We then 
discuss how these values of t> lt and D' %t lead to a wrong 
estimate of the contact polarization before comparing our 
results for r^ lt and D^ lt qualitatively with measurements 
on graphene in the literature. The modeled effect of the 
contacts on spin transport only depends on the resistance 
of the barrier and not on the type of barrier. Hence, al- 
though most contact interfaces used in the non-local ge- 
ometry to study spin transport in graphene are not truly 
in the tunneling regime, we can conclude that with the 
the resistance of the commonly used barriers the effect of 
back diffusion into the contacts on the spin transport is 
only minor and the spi n transport properties a re mainly 
limited by other effects ! 5 ' 7 ' 10 ' 14 ' 1 ™ 20 ' 33 ' 37 ' 39 ' While ex- 
plicitly discussing the effect of low resistive contacts on 
the non-local geometry, similar effects also play a role for 
local measurements^ 7 ? 

We also briefly discussed the guidelines for a good and 
reliable Hanle fit as an incorrectly performed fit can also 
lead to misinterpretations of the spin transport proper- 
ties of a diffusive channel while a correct fit leads to un- 
ambiguous results. 
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The discussion in the main text focused on the sym- 
metric case when the injector and the detector contacts 
have equal i?-parameters. Here we address the gen- 
eral case of dissimilar injector and detector contacts and 
demonstrate an equivalence that allows us to map this 
general case to the more symmetric one presented above. 

The general expression for the non-local resistance R n i 
given bjSsl 



Rni = 



RsnX 



2 / p 



2R t 



_p2 



n h+i 



x 



P 2 



- exp (-2L/X) 



.1 (Al) 



where R\^ correspond to the i?-parameters of the injec- 
tor and detector contacts, and the rest of the parame- 
ters are the same as those presented in the discussion of 

E q . (fur 

This equation can be simplified by realizing that for 
highly spin polarized contacts (P ~ 1) there is no contact 
induced spin relaxation, even for low resistance contacts. 
Therefore if we consider 1 — P 2 ss 1 we obtain 
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P 2 RsaX 



2R 1 
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exp (—L/X) 



2W 
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2# 2 
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which has a similar structure as Eq. ([l]) 19 . Following 
simple algebra, we can equate both equations and solve 
for the i?-paramctcr of Eq. ([I]) , which can be understood 
as an effective i?-parameter R e ff(Ri, R 2 , L, A) given by, 
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2 It; 
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FIG. 6: (Color online) Mapping the problem of dissimilar 
contacts Ri ^ Ki into the simpler one of identical contacts 
with a common R e ff- (a) Two-dimensional map of equiva- 
lent R e ff/X as a function of Ri/X of the contacts for L/X = 1 
from Eq. (A3 1, (b) Normalized deviation of i? e ///A obta ined 
from Eq. (A4| relative to the exact result from Eq. (A3 1 
in the limit of L/X 3> 1. The values are normalized using 
(R e f f (L/X = 10) - R eff (L/X = 0))/(R eff (L/X = 10). 



allowing us to map the case of dissimilar contacts into the 
symmetric case of equal contacts with — Reff - One 
example of such a mapping is shown in Fig. [6] (a) for the 
representative case of L/X — 1 and Ri/X = 0.1-10. We 
observe that when R\ ^ R2 then R e ff ~ min (Ri, R2) 
and also the trivial case of R\ = R2 = R e f / ■ 

The exact mapping depends on L and on A, which 
requires careful application to analyze experimental data. 
We remark th at th is issue is absent for the case of 2L/A sw 
0, where Eq. (A3) reduces to the simple form, 



R 



eff 



2 \Ri 



1 

R~2 



(A4) 



equivalent to a 3-terminal measurement where both con- 
tacts are in a parallel configuration. 



Although Eq. ( A4 ) is strictly speaking valid only when 



both contacts are closely spaced, we have observed that 
it offers a reasonable approximation even at finite sepa- 
ration L between the contacts. In Fig. [6] (b) we compare 
the resulting R e ff/X for the extr eme case of large sepa- 
ration (2L/X ^> 1) from Eq. (A3), to the value obtained 
from the simpler Eq. (A4). Surprisingly, in the experi- 



mentally relevant range of intermediate conductivity mis- 
match Ri/X 



0.3-10, Eq. (A4 



result at infinite separation on 
a strong conductivity mismatch (Ri/X < 0.1 
apply the exact result of Eq. (A3). 



) deviates from the exact 
y by less than 20%. For 
one should 



Appendix B: Fitting simulated Hanle curves 

The research presented in this paper is based on the 
following concept: Hanle precession curves are simulated 



B(T) 
50 100 150 



BCD 




R/A=0.01, L/A=0.1 

sim. data 

- fit 

R/A=0.1, L/A=3 

sim. data 

fit 



0.1 0.2 0.3 
B(T) 




B(T) 



FIG. 7: (Color online) Two sets of simulated data with the 
corresponding fits assuming an amplitude of 1 and a baseline 
at 0. Both sets were simulated for D = 0.01 m 2 /s, r = 100 ps 
and the values for L/X and R/X shown in the legend, (a) 
shows the full Hanle curve, (b) zooms in on the part close to 
B = 0. The _B-scale for the curve with the red dashed fit is on 
the bottom, the scale for the curve with the light blue dotted 
fit on top. 



following the model presented in Ref. [19] including con- 
tact induced spin relaxation and are fitted neglecting the 
contact induced effects. Fig. [7] shows how well the simu- 
lated data can be fitted with a Hanle curve for different 
values of R/X and L/X. The curve for R/X = 0.1 and 
L/X = 3 shows that even for small R/X- values (corre- 
sponding to a contact resistance of R = 100 when 
R sq — 1 kf2 and W = 1 fim) we get an excellent fit 
(although with a reduced r^ lt and incr eased D* u ). On 
the other hand the curve and fit for the combination 
R/X = 0.01 and L/X = 0.1 points out that the fit is 
not describing the curve properly for very small values of 
the two parameters. This is especially well visible close 
to B = (Fig. [7] (b)) where due to the strong contact in- 
duced relaxation we observe a distinct drop of R n i which 
the fit cannot describe. In Fig. [3] (c) and (d) of the main 
text it is visible for which sets of parameters the fits do 
not describe the Hanle curves well, as those are the points 
that do not show a smooth line shape when plotting D^ lt 
MS <X function of R/X or L/X. 



Appendix C: r /lt in the limit L/X < 1 

We can obtain T^ lt by performing a standard Hanle fit 
on simulated data that includes contact induced relax- 
ation. Here we show that we can approximate the value 
of T^ %t for small L using an easy reasoning. 
In the limit L/A< 1 our system resembles the 3-terminal 
Hanle geometry^ as we have two contacts connected to 
approximately the same point of the transport channel 
with one of the contacts injecting spins and the other 
detecting them. At the same time there is the trans- 
port channel pointing in two directions away from the 
injection point. Therefore we get for the spin resistance 
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1 / R*s P in = 1 / R s c h + 1 / Rc ■ K we now take the ratio of the 
spin resistance including the contact resistance (R* s 



and R S pi n — R ch 



(Kpin for R c -> °°) we get: 



R 



spin 



l/R 



S p 'i l I 



R M 



l/R 



Spll I 



_ R/X 

1 / Rs ch + l l R c ~ 1 + 



spin ) 
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The spin resistance is proportional to the non-local signal 
(Rspin oc R n i) and for L = and R — > oo the non- 
local signal is proportional to the spin relaxation length 
R n i oc A (see Eq. ([!])). Therefore we get: 



R* spin oc \f lt we have to keep in mind that X^ lt is ob- 
tained assuming R — > oo so we have to assume this also 
here in this analysis. 

The relation between the ratio of the spin relaxation 
lengths and the ratio of the spin relaxation times is valid 
as D^/D « 1. 
Hence, we get the result: 



-fit 



R/X 
1 + R/X 



(C3) 



spin 

Rs 13171 



\fit 



T fit 
T 



(C2) 



We can use here for both R* pin and R sp i n the relation 
R n i oc A. This is obviously valid for R sp i n and for 



In the limit L/X <C 1 we therefore expect 
T f lt /T = (0.98,0.83,0.25,8.3 x 10- 3 ,9.8 x 10~ 5 ) 
for R/X = (100,10,1,0.1,0.01) in good agreement with 
the values in Fig. [3] (b) in the limit L/X <C 1. 



* Electronic address: t.maassen@rug.nl 

1 S. A. Wolf, A. Y. Chtchelkanova, and D. M. Treger, IBM 
J. Res. Dev. 50, 101 (2006). 

2 J. Fabian, A. Matos-Abiague, C. Ertler, P. Stano, and 
I. Zutic, Acta Phys. Slov. 57, 565 (2007). 

3 G. Schmidt, D. Ferrand, L. W. Molenkamp, A. T. Filip, 
and B. J. van Wees, Phys. Rev. B 62, R4790 (2000). 

4 E. I. Rashba, Phys. Rev. B 62, R16267 (2000). 

5 N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and 
B. J. van Wees, Nature 448, 571 (2007). 

6 D. Huertas-Hernando, F. Guinea, and A. Brataas, Eur. 
Phys. J. Spec. Top. 148, 177 (2007). 

7 C. Jozsa, T. Maassen, M. Popinciuc, P. J. Zomer, 
A. Veligura, H. T. Jonkman, and B. J. van Wees, Phys. 
Rev. B 80, 241403 (2009). 

8 C. Ertler, S. Konschuh, M. Gmitra, and J. Fabian, Phys. 
Rev. B 80, 041405 (2009). 

9 A. H. Castro Neto and F. Guinea, Phys. Rev. Lett. 103, 
026804 (2009). 

T.-Y. Yang, J. Balakrishnan, F. Volmer, A. Avsar, 
M. Jaiswal, J. Samm, S. R. Ali, A. Pachoud, M. Zeng, 
M. Popinciuc, et al., Phys. Rev. Lett. 107, 047206 (2011). 

1 P. Zhang and M. W. Wu, New J. Phys. 14, 033015 (2012). 

2 H. Ochoa, A. H. Castro Neto, and F. Guinea, Phys. Rev. 
Lett. 108, 206808 (2012). 

3 K. Pi, W. Han, K. M. McCreary, A. G. Swartz, Y. Li, and 
R. K. Kawakami, Phys. Rev. Lett. 104, 187201 (2010). 

4 T. Maassen, F. K. Dejene, M. H. D. Guimaraes, C. Jozsa, 
and B. J. van Wees, Phys. Rev. B 83, 115410 (2011). 

5 T. Maassen, J. J. van den Berg, N. IJbema, F. Fromm, 
T. Seyller, R. Yakimova, and B. J. van Wees, Nano Lett. 
12, 1498 (2012). 

6 M. H. D. Guimaraes, A. Veligura, P. J. Zomer, T. Maassen, 
I. J. Vera-Marun, N. Tombros, and B. J. van Wees, Nano 
Lett. 12, 3512 (2012). 

7 K. M. McCreary, A. G. Swartz, W. Han, J. Fabian, and 
R. K. Kawakami, arXiv:1206.2628v2 (2012), 1206.2628. 

8 P. J. Zomer, M. H. D. Guimaraes, N. Tombros, and B. J. 
van Wees, arXiv:1209.1999vl (2012). 

9 M. Popinciuc, C. Jozsa, P. J. Zomer, N. Tombros, 
A. Veligura, H. T. Jonkman, and B. J. van Wees, Phys. 



Rev. B 80, 214427 (2009). 

W. Han, K. Pi, K. M. McCreary, Y. Li, J. J. I. Wong, 

A. G. Swartz, and R. K. Kawakami, Phys. Rev. Lett. 105, 
167202 (2010). 

B. Dlubak, P. Seneor, A. Anane, C. Barraud, C. Deranlot, 
D. Deneuve, B. Servet, R. Mattana, F. Petroff, and A. Fert, 
Appl. Phys. Lett. 97, 092502 (2010). 

B. Dlubak, M.-B. Martin, C. Deranlot, B. Servet, 
S. Xavier, R. Mattana, M. Sprinkle, C. Berger, W. A. 
De Heer, F. Petroff, et al., Nat. Phys. 8, 557 (2012). 
J. Abel, A. Matsubayashi, T. Murray, C. Dimitrakopoulos, 
D. B. Farmer, A. Afzali, A. Grill, C. Y. Sung, and V. P. 
LaBella, J. Vac. Sci. Technol. B 30, 04E109 (2012). 
F. J. Jedema, H. B. Heersche, A. T. Filip, J. J. A. Basel- 
mans, and B. J. van Wees, Nature 416, 713 (2002). 
N. Tombros, S. J. van der Molen, and B. J. van Wees, 
Phys. Rev. B 73, 233403 (2006). 

S. Takahashi and S. Maekawa, Phys. Rev. B 67, 052409 
(2003). 

For a negligible spin resistance of the ferromagnetic elec- 
trodes compared to the one of the barrier and to the resis- 
tance of the transport channel (R sq ) we call the resistance 
of the interface barrier between the contact electrode and 
the transport channel "contact resistance" (Rc) - The con- 
tact resistance consists of the parallel resistance for spin 
up (i?t) and spin down (iij.). We assume Rf as R± as we 
are discussing polarizations of P < 30% leading to an er- 
ror in determining the effect of the contact resistance on 
the spin relaxation of less than 10%. In the case of trans- 
parent contacts (i.e. no high resistive barrier), Rc is equal 
to the spin resistance of the ferromagnetic leads Rf pin . In 
this case we also have to replace the spin polarization of 
the barrier P by the spin polarization of the ferromagnetic 
leads P F . 

We are discussing here a 2-dimensional channel like 
graphene. In case of a 3-dimensional channel R would be 
R — (Rc / p)A where A is the cross section and p the resis- 
tivity of the transport channel. 

M. Johnson and R. H. Silsbee, Phys. Rev. B 37, 5326 
(1988). 

F. J. Jedema, M. V. Costache, H. B. Heersche, J. J. A. 



9 



Baselmans, and B. J. van Wees, Appl. Phys. Lett. 81, 5162 
(2002). 

To fit the data we first normalize the simulated R n i values 
with R n i(B = 0) and keep the fitting amplitude fixed as 
1. This enables a more standardized fitting procedure and 
reduces the fitting time considerably. The fixed amplitude 
can be justified as fitting with a non-fixed amplitude gives 
in most cases a variation of the amplitude by less than 1% 
and only varies more in the regime where R/X < 0.1 and 
L/X < 0.1. 

S. P. Dash, S. Sharma, R. S. Patel, M. P. de Jong, and 
R. Jansen, Nature 462, 491 (2009). 

M. Wojtaszek, I. J. Vera-Marun, T. Maassen, and B. J. 
van Wees, arXiv:1209.2365vl (2012). 

That the data agrees only qualitatively with our simula- 
tions and that the second data set does not show a clear 
dependence of r on L can be related to differences in the 
i?-parameter between the electrodes as Wojtaszek et al. 



report ^-parameters of 3 /im < R < 100 fim. Electrically 
floating contacts between the injector and detector elec- 
trode could also be influencing the data. 
T. Suzuki, T. Sasaki, T. Oikawa, M. Shiraishi, Y. Suzuki, 
and K. Noguchi, Appl. Phys. Express 4, 023003 (2011). 
N. Tombros, S. Tanabe, A. Veligura, C. Jozsa, M. Popin- 
ciuc, H. T. Jonkman, and B. J. van Wees, Phys. Rev. Lett. 
101, 046601 (2008). 

I. J. Vera-Marun, V. Ranjan, and B. J. van Wees, Nat. 
Phys. 8, 313 (2012). 

T. Maassen, J. J. van den Berg, E. H. Huisman, H. Di- 
jkstra, F. Fromm, T. Seyller, and B. J. van Wees, 
arXiv:1208.3129vl (2012). 

W. Han and R. K. Kawakami, Phys. Rev. Lett. 107, 
047207 (2011). 

H. Jaffres, J.-M. George, and A. Fcrt, Phys. Rev. B 82, 
140408 (2010). 



